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Abstract 



We have measured the cosmic ray spectrum at energies above 10 eV using the two 
air fluorescence detectors of the High Resolution Fly's Eye experiment operating 
in monocular mode. We describe the detector, PMT and atmospheric calibrations, 
and the analysis techniques for the two detectors. We fit the spectrum to models 
describing galactic and extragalactic sources. Our measured spectrum gives an ob- 
servation of a feature known as the "ankle" near 3 x 10 18 eV, and strong evidence 
for a suppression near 6 x 10 19 eV. 



1 Introduction 

The highest energy cosmic rays yet 
detected, of energies up to and above 
10 20 eV, are interesting in that they 
shed light on two important ques- 
tions: how are cosmic rays acceler- 
ated in astrophysical sources, and 
how do they propagate to us through 
the cosmic microwave background 
radiation (CMBR)[1]? The acceler- 
ation of cosmic rays to ultra high 
energies is thought to occur in exten- 
sive regions of high magnetic fields, 
regions which are expanding at rela- 
tivistic velocities [2]. Such regions are 
rare and are to be counted among the 
most violent and interesting objects 
in the universe. 

Once accelerated, interactions be- 
tween the ultra high energy cos- 
mic rays (UHECR) and the CMBR 
cause the cosmic rays to lose energy. 
The strongest energy loss mecha- 
nism comes from the production of 
pions in these CMBR interactions 
at UHECR energies above about 
6 x 10 19 eV. This energy loss mecha- 
nism produces the Greisen-Zatsepin- 
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K'uzmin (GZK) suppression [3, 4]. In 
addition, e + e~ production in these 
same interactions provides a some- 
what weaker energy loss mechanism 
above a threshold of about 5 x 10 17 
eV. A third important energy-loss 
mechanism at all energies comes 
from universal expansion. 

In previous publications[5,6] , we have 
reported on our measurements of the 
cosmic ray spectrum using data col- 
lected independently, in monocular 
mode, by the two detectors of the 
High Resolution Fly's Eye experi- 
ment (HiRes). We here report on an 
updated measurement of the flux of 
UHECR, covering an energy range 
from 2.5 x 10 17 eV to over 10 20 eV, 
using a significantly larger data set 
for the HiRes-II detector. With the 
improved statistical power available 
in this data, we study two features in 
this spectrum: a break in the spec- 
tral slope at 3 x 10 18 eV, called the 
"ankle" [7,8,9], and a steepening of 
the spectrum near the threshold for 
pion production. 

The HiRes experiment performs a 
calorimetric measurement of the en- 
ergy of cosmic rays. UHECR pro- 
duce extensive air showers (EAS) 
when they enter the atmosphere. The 
HiRes detector collects the fluores- 
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cence light emitted by EAS as they 
propagate through the atmosphere. 
Charged particles in the shower ex- 
cite nitrogen molecules which fluo- 
resce in the ultraviolet (300 to 400 
nm). The fluorescence yield is about 
five photons per minimum ionizing 
particle per meter of path length [10]. 
As an EAS propagates through the 
atmosphere, the detector measures 
the number of photons seen as a 
function of time and angle. From 
this information, we reconstruct the 
geometry of the shower and the solid 
angle subtended by the detector 
from each point of the shower. From 
the number of photons collected, we 
reconstruct the number of charged 
particles in the shower as a function 
of the depth of the atmosphere tra- 
versed. We integrate the energy de- 
posited in the atmosphere [11] to find 
the energy of the primary cosmic ray. 



UHECRs are thought to be protons 
or heavier nuclei up to iron. While 
nucleus-nucleus collisions are com- 
plex, the general features of the in- 
teraction can be understood in terms 
of a simple superposition model. In 
this model each nucleon generates 
an independent EAS. The superpo- 
sition of many, lower energy showers 
will result in an EAS with different 
statistical properties than an EAS 
produced by one high energy proton. 
This allows one to measure the com- 
position of the primary cosmic rays 
on a statistical basis. 



2 The HiRes Detectors 



The HiRes detectors have been de- 
scribed extensively elsewhere [12, 13]. 
In brief, they consist of spherical 
mirrors, of area 5.1 m 2 , which collect 
the fluorescence light and focus it 
onto a cluster of 256 photomultiplier 
tubes arranged in a 16 x 16 array. 
Each tube in the cluster views about 
one square degree of the sky. Time 
and pulse height information are col- 
lected from each tube. The HiRes 
detectors trigger on and reconstruct 
showers that occur within a radius 
of about 35 km. 

The HiRes-I detector is located atop 
Little Granite Mountain on the U.S. 
Army Dugway Proving Ground in 
west-central Utah. It consists of 21 
mirrors, and their associated photo- 
tube arrays, arranged in one ring, 
observing from 3 degrees to 17 de- 
grees in elevation and providing al- 
most complete coverage in azimuthal 
angle. The detector uses a sample- 
and-hold readout system which inte- 
grates phototube pulses for 5.6 fis. 
This is long enough to collect the 
signal from all cosmic ray showers of 
interest. 

The HiRes-II detector is located on 
Camel's Back Ridge, also on Dug- 
way Proving Ground, about 12.6 km 
SW of HiRes-I. It consists of 42 mir- 
rors, arranged in two rings, covering 
from 3 to 31 degrees in elevation 
and almost the whole azimuthal an- 
gle range. This detector uses a flash 
ADC (FADC) readout system with 
a 100 ns sampling time. 
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In this article, we present data col- 
lected from June 1997 to February 
2003 for HiRes-I, and from Decem- 
ber 1999 through September 2001 
for HiRes-II. For HiRes-II, this is 
about four times the data that was 
reported on previously. We collect 
data on nights when the moon is 
down for three hours or more. In a 
typical year each detector collects up 
to about 1000 hours of data. 

The weather is clear about 2/3 of 
the time at the HiRes sites. Since 
clouds can reduce the experiment's 
aperture, we record the existence of 
clouds by operator observations, in- 
frared cameras, and evidence from 
data collected by the detector (this 
consists primarily of the upward go- 
ing laser and flasher pulses, used to 
measure the atmospheric conditions, 
which have a distinct signature upon 
encountering a cloud; actual cos- 
mic rays also appear emerging from 
clouds). Only data from those nights 
in which the aperture is not reduced 
by cloud cover are used in our spec- 
trum measurements. 



3 Calibration 

The two most important calibrations 
we perform are of the photomulti- 
plier tube (PMT) gains[14,15], and 
of the clarity of the atmosphere [16]. 
We use a stable xenon flash lamp, 
carried to each detector and used to 
illuminate the photomultiplier array, 
to find PMT gains. The xenon lamp 
produces a light intensity of about 
10 photons per mm 2 at the face of 
the PMTs; this intensity is traceable 



to NIST-calibrated photodiodes and 
is stable to about 2%, flash-to-flash. 
Separate calibrations of PMT gains 
using photoelectron statistics and 
using the absolute light intensity of 
the xenon flash lamp agree within 
uncertainties. Xenon flash lamp data 
are collected about once a month. 
A second calibration system, using 
a frequency-tripled YAG laser, is 
used to monitor phototube gains on 
a night-to-night basis. We estimate 
that the relative calibration tech- 
niques are accurate to about 3% with 
an absolute calibration uncertainty 
of about ±10%. 

The atmosphere is our calorimeter, 
but it is also the medium through 
which fluorescence light propagates 
to the detectors. To calculate the 
number of fluorescence photons emit- 
ted by a cosmic ray shower, we must 
understand the way in which the 
atmosphere scatters this light be- 
tween the EAS and the detector. The 
molecular component of the atmo- 
sphere is quite constant, with only 
small seasonal variations, and the 
Rayleigh scattering it produces is 
well understood. The aerosol content 
of the atmosphere can vary consid- 
erably over time, and with it, the 
amount of light scattered and its 
angular distribution. 

To measure these quantities, we per- 
form an atmospheric calibration us- 
ing YAG lasers operating at wave- 
length A = 355 nm. At each of our 
two sites, we have a steerable beam 
laser which is fired in a pattern of 
shots that covers the detector's aper- 
ture, and which is repeated every 
hour. The scattered light from the 
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laser at one site is collected by the de- 
tector at the other site. The amount 
of detected light is then analyzed to 
determine the scattering properties 
of the atmosphere. The properties 
that we measure are the vertical 
aerosol optical depth (VAOD), the 
horizontal aerosol extinction length, 
and the aerosol scattering phase 
function (the angular distribution of 
the differential scattering cross sec- 
tion). 

Because about half of the data from 
HiRes-I were collected before the 
lasers were installed, we use aver- 
age values of the measured param- 
eters in this analysis: a horizontal 
aerosol extinction length of 25 km 
(the average horizontal molecular 
extinction length is 17 km), an av- 
erage phase function, and a VAOD 
of 0.04[5,16,17]. The atmosphere at 
our sites is quite clear: the average 
atmospheric correction to an event's 
energy is about 10% (see below for 
the effect on flux measurements). 
We are most sensitive to the value 
of VAOD. The RMS of the VAOD 
distribution is 0.02, and we use this 
RMS value conservative esti- 

mate of the systematic uncertainty 
in this parameter. 



4 HiRes-II Data Analysis 

The analysis of the HiRes-II monocu- 
lar data has been described previously[6] . 
The data presented here were col- 
lected during 540 hours of good 
weather running, and consists of 21 
million triggers, mostly of random 
sky noise and events generated by 



atmospheric lasers and other man- 
made light sources. Events were se- 
lected that satisfied the following 
criteria: 

• Angular speed < 11° /is -1 

• Selected tubes > 6 

• Photoelectrons/degree > 25 

• Track length > 7°, or > 10° 

for events extending above 
17° elevation 

• Zenith angle < 80° 

• In-plane angle < 130° 

• In-plane angle uncertainty < 

30° 

• 150 < X max < 1200 g/cm 2 , 

and is within 50 g/cm 2 of 
begin visible in detector 

• Average Cerenkov Correc- 

tion < 70% 

• Geometry fit x 2 /d.o.f. < 10 

• Profile fit x7d-o.f. < 10 

• Minimal trigger from signal 

tubes required after March 
2001 

These cuts remove events in which 
the monocular geometric reconstruc- 
tion is poor or in which the longitu- 
dinal profile cannot be determined 
accurately. The final event sample 
consisted of 2685 events covering 
an energy range from 1.6 x 10 17 eV 
(log 10 E = 17.2) to 10 20 eV. 

The geometry of each event is re- 
constructed using the time and an- 
gle information from the hit PMTs. 
First a pattern recognition step is 
performed to choose phototubes that 
lie on a line both in angle and in 
time. Next the plane that contains 
both the shower and the detector is 
determined from the azimuth and el- 
evation of hit tubes; the angle of the 
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shower in this plane is determined 
from a fit to phototube time and 
angle information. The resolution of 
shower-detector plane determination 
is about 0.6°, and the in-plane angle 
uncertainty is 5° on average. 

With the geometry determined, the 
profile of the number of charged par- 
ticles in the shower is calculated from 
the phototube pulse heights. Cor- 
rections are made for atmospheric 
scattering of the light, and for other 
effects such as mirror reflectivity, 
phototube quantum efficiency, etc. A 
correction is made for the Cerenkov 
light produced by charged particles 
in the shower. Both direct and scat- 
tered Cerenkov light contributions to 
the light seen by the PMTs are calcu- 
lated and subtracted. The number of 
charged particles is calculated from 
the fluorescence light at the shower 
using the fluorescence yield and its 
pressure and temperature variation 
as given by Kakimoto et al[10]. The 
resulting shower development pro- 
file, expressed as a function of slant 
depth, is fit to the Gaisser-Hillas 
parameterization [18] (this has been 
seen to fit UHE cosmic ray showers 
quite well [11, 19]). We integrate over 
the Gaisser-Hillas function and mul- 
tiply by the average energy loss rate 
of 2.19 MeV/g/cm 2 to calculate the 
energy of the primary cosmic ray. We 
then correct for unobserved energy, 
mostly neutrinos and muons which 
hit the ground. This correction[17], 
which varies from 10% at 3 x 10 17 
eV to 5% at 10 19 eV, is determined 
during the Monte Carlo calculation 
of the aperture. It is similar to the 
calculation in reference [11]. 
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Fig. 1. A scatter plot of the HiRes-II 
energy calculated using monocular ge- 
ometry versus the energy calculated us- 
ing the stereo geometry for those events 
observed in stereo. 

A fraction of the HiRes-II events 
are also observed by HiRes-I. In this 
case we perform a cross-check on 
our monocular determination of the 
shower geometry. Figure 1 shows 
a scatter plot of the energy using 
monocular geometry versus the en- 
ergy using stereo geometry, for those 
events in which such a comparison is 
possible. 

The energy resolution, including sta- 
tistical and systematic effects and 
the uncertainty in reconstructing the 
shower geometry, has been calculated 
in the Monte Carlo simulation [17]. 
The overall resolution is ±17%. It 
improves from ±18% below 10 18 eV 
to ±12% above 10 19 eV. 



5 Monte Carlo Simulation 

To calculate the aperture as a func- 
tion of cosmic ray energy, a very ac- 
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curate Monte Carlo (MC) simulation 
of the experiment was performed [6]. 
Two libraries of cosmic ray showers, 
one for proton primaries and one for 
iron primaries, were generated using 
the Corsika 5.61 [20] EAS simula- 
tion program and the QGSJet 01 [21] 
hadronic event generator. Events 
from these libraries were placed by 
a detector simulation in the vicinity 
of the HiRes-II detector. This pro- 
gram also simulated the fluorescence 
and Cerenkov light generated by the 
showers, and calculated how much of 
this light would have been collected 
by the detectors. A complete simula- 
tion of the optical path, trigger, and 
readout electronics was performed. 
This simulation followed the exper- 
imental conditions that pertained 
over the data- collect ion period. The 
results were written out in the same 
format as the data and analyzed 
by the same data analysis program 
described above. The stereoscopic 
energy spectrum of the Fly's Eye 
experiment [7], in the form of a bro- 
ken power law fit, and the compo- 
sition measurements made by the 
HiRes/MIA hybrid experiment [22] 
and by HiRes in stereo [23] were used 
as inputs. 

To convince ourselves that the MC 
simulation is accurate, we compare 
many MC distributions of geomet- 
rical and kinematic variables to the 
data. The agreement in these com- 
parisons is excellent and indicates 
that we understand our detector. Fig- 
ure 2 shows the brightness of show- 
ers: the number of photoelectrons 
per degree of track. The agreement 
between the data and MC simulation 
indicates that the same amount of 




Fig. 2. Comparison of HiRes-II data 
and MC for the photoelectrons per de- 
gree of track. In the upper frame, the 
filled squares with the histogram are the 
data, the open squares are the MC. The 
lower frame shows the ratio of data to 
MC for each bin. 

light is collected in the MC as in the 
data. Figure 3 shows the \ 2 °f a fit 
to the time vs. angle plot from which 
we determine shower geometry. The 
agreement here indicates that the 
resolution of the MC is the same as 
that of the data. Figure 4 shows a 
histogram of the number of events 
vs the logarithm of their energy in 
EeV. The agreement here shows that, 
when we use previous measurements 
of the spectrum and composition in 
the MC, and a complete simulation 
of the acceptance, we reproduce the 
experiment's energy dependence. 



6 HiRes-I Analysis 



The analysis of the HiRes-I monoc- 
ular data has also been described 
previously [5,6]. The main difference 
from the HiRes-II analysis is that, 
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X 2 / DOF 



Fig. 3. Comparison of HiRes-II data and 
MC for the x 2 of a fit to the time vs an- 
gle plot assuming a vertical shower. In 
the upper frame, the filled squares with 
the histogram are the data, the open 
squares are the MC. The lower frame 
shows the ratio of data to MC for each 
bin. 
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Fig. 4. Comparison of HiRes-II data 
and MC for the reconstructed energy. In 
the upper frame, the filled squares with 
the histogram are the data, the open 
squares are the MC. The lower frame 
shows the ratio of data to MC for each 
bin. 



with only one ring of mirrors, most 
tracks are too short to reliably de- 
termine the geometry from timing 
alone. Although the determination 
of the shower-detector plane is still 
excellent, correlations between the 
fit distance to the shower and the fit 
in-plane angle become large for short 
tracks. 

A reconstruction procedure using 
the pulse height information in ad- 
dition to the tube angles and timing 
information has been developed: the 
profile constrained fit (PCF). The 
PCF uses the one-to-one correlation 
between in-plane angle and shower 
profile: the in-plane angle with the 
best fit shower profile is chosen as 
the in-plane angle of the shower. The 
Gaisser-Hillas function is used in this 
profile fit. The PCF works poorly for 
events close to the detector (within 
about 5 km), and for lower energy 
events (below 3 x 10 18 eV), where less 
of the shower profile is seen. These 
events are excluded from the HiRes- 
I monocular sample. The PCF also 
works poorly if too much Cerenkov 
light contaminates the fluorescence 
signal; these events are cut also. In 
reconstructing MC events, it is found 
that, even with these cuts, the res- 
olution is somewhat worse than for 
HiRes-II, and that there is an energy 
bias. 

Since stereo events are seen in both 
detectors, they have excellent geo- 
metrical determination using the in- 
tersection of the two shower-detector 
planes. For these events, compari- 
son of the PCF reconstruction to the 
stereo reconstruction shows the same 
energy resolution and bias as seen in 
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Fig. 5. Energy resolution using PCF, 
after bias correction. The histogram 
shows MC resolution, the data points 
show the data monocular resolution in 
stereo events. For the MC, Erec refers 
to the reconstructed, monocular energy, 
while Emc refers to the generated en- 
ergy in the same events. For the data, 
Emono refers to the energy recon- 
structed using the monocular geometry 
(corresponding to Erec m the MC), 
while Est refers to the energy recon- 
structed using the stereo geometry. 

the MC sample. Having confidence 
that we understand the PCF, we 
correct for the bias. Figure 5 shows 
the energy resolution of the PCF re- 
construction for MC events and for 
stereo events after the correction. 
The agreement is excellent. 

Figure 6 shows comparisons between 
the HiRes-I data and the MC simu- 
lation for the distance to the shower 
core of showers in three energy bands. 
Again the agreement is excellent. 



7 Systematic Uncertainties 



The largest systematic uncertainties 
in the calculations of energy are the 
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Fig. 6. Shower core distance distribu- 
tion using PCF, in HiRes-I data and in 
MC. The squares with error bars are the 
data, the histogram is the MC. 

absolute calibration of the photo- 
tubes (±10%) [15], the fluorescence 
yield (±10%) [10], and the correction 
for unobserved energy in the shower 
(±5%) [11,17]. These three uncer- 
tainties, added in quadrature, give an 
uncertainty in the energy of ±15%. 
This effect of this energy uncertainty 
in calculating the flux is ±27%[5]. 

To test the sensitivity of the flux mea- 
surement to atmospheric uncertain- 
ties, we generated new MC samples 
with VAOD values of 0.02 and 0.06, 
i.e., with the average plus and minus 
one RMS value, and analyzed them 
(and the data) using the same VAOD 
values. This provides a conservative 
estimate of the flux uncertainty since 
the systematic uncertainty in the av- 
erage VAOD is less than the RMS. 
The result was a change in the flux 
of ±15%. Adding this in quadrature 
with the sources of systematic uncer- 
tainty described above results in a 
net uncertainty of ±31%. This uncer- 
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tainty is common to the flux measure- 
ments from HiRes-I and HiRes-II. 

The effect of using an average VAOD 
value, rather than the changing 
but measured values has also been 
studied [24, 17]. The changes to the 
energy scale and flux are negligible. 

The limited elevation coverage of the 
HiRes-II detector makes the aperture 
calculation sensitive to the compo- 
sition assumptions used in the MC 
simulation. This and other sources of 
systematic uncertainty (the given in- 
put spectrum and using an average 
atmosphere) are considered in refer- 
ence [24]. The composition assump- 
tions have a negligible effect on the 
aperture above an energy of 10 18 eV, 
and give a systematic uncertainty of 
order the statistical uncertainty only 
at 3 x 10 17 eV. 



8 Results 

Figure 7 shows the calculated aper- 
ture of the two HiRes detectors. At 
an energy of 10 20 eV the aperture is 
nearly 10,000 km 2 sr. 

Figure 8 shows the measured spec- 
trum of cosmic rays[25]. The spec- 
trum has been multiplied by E 3 for 
clarity. The closed squares (open cir- 
cles) are the HiRes-I (HiRes-II) mea- 
surements. For comparison to pre- 
vious experiments, the up-triangles 
are the stereo Fly's Eye spectrum[7], 
and the down-triangles are the result 
of the Akeno Giant Air Shower Ar- 
ray (AG AS A) [8]. The HiRes-I and 
HiRes-II monocular measurements 
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Fig. 7. The calculated apertures of 
HiRes-I (red squares) and HiRes-II 
(black circles) as a function of energy. 

agree with each other very well in the 
overlap region, and are also in good 
agreement with the Fly's Eye stereo 
spectrum. 

In this plot the ankle shows up clearly 
at 3 x 10 18 eV (log 10 E = 18.5). The 
spectrum steepens again at 6 x 10 19 
eV (log 10 E = 19.8). The AGASA 
spectrum appears to continue un- 
abated above this energy while the 
HiRes spectrum falls above this 
point. 

We test whether our data are consis- 
tent with this interpretation of the 
AGASA spectrum by fitting our data 
to a broken power law. This fit is 
also shown on Figure 8. This fit had 
two floating break points separating 
three regions of constant spectral 
slope. The fit was performed using 
the normalized, binned maximum 
likelihood method[26], which allows 
us to include sparsely populated and 
empty bins. The fitted break points 
are at log 10 E (in eV) = 18.47 ± 0.06 
and 19.79 ± 0.09. The fitted spec- 
tral slopes are 7 = 3.32 ± 0.04, 
2.86 ± 0.04, and 5.2 ± 1.3. The \ 2 
for the fit is 30.1 for 33 degrees of 
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Fig. 8. -E 3 times the UHECR Flux. Results from the HiRes-I (red squares) and 
HiRes-II (black circles) detectors, the AGASA experiment (blue down-triangles) and 
the Fly's Eye experiment (in stereo mode; magenta up-triangles) are shown. The 
Fly's Eye points have been shifted to the right by one quarter bin (A log 10 E = 0.025) 
for clarity. Also shown are two spectral law fits to the HiRes-I and HiRes-II spectra 
as described in the text. The la upper limits for two empty bins of each HiRes 
spectra are also shown. 



freedom. If we extend the middle 
section of the fit (as shown by the 
red/grey line in Figure 8) to higher 
energies, our aperture predicts that 
we should have 28.0 events above 
the second break point energy of 
log 10 -E = 19.79, where we really 
have 11. The Poisson probability for 
11 or fewer events with a mean of 
28 is 2.4 x 10" 4 . We therefore con- 
clude that our data is not consistent 
with a continuation of the spectrum 
unabated above the pion production 
threshold. It is worth emphasizing 
that we have considerable sensitivity 
to such a continuation, but the data 
do not support it. 



A similar fit with only one break 
point has a x 2 °f 46.0 for 35 de- 
grees of freedom, worse by nearly 16 
than the fit above (a significance of 
~ 3.7er). The fitted break point is at 
log 10 E = 18.45 ±0.03, and the fitted 
spectral slopes are 7 = 3.32 ± 0.03 
and 2.85 ± 0.05. A fit with no break 
points at all has a bad \ 2 °f H4 for 
37 degrees of freedom, demonstrat- 
ing that the ankle is clearly observed 
in our data. The spectral slope is 
7 = 3.12 ± 0.01. This fit is shown on 
Figure 8 as a cyan/light grey line. 
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9 Fitting the Spectrum 

The implications of our spectrum 
measurement can be explored using a 
toy model of UHECR. In this model, 
there are two types of sources, galac- 
tic and extragalactic. We choose the 
galactic sources to be consistent with 
the HiRes/MIA and HiRes stereo 
composition measurements [22, 23]: 
we assign the iron component of the 
cosmic ray flux to be galactic [27]. 
This assignment is consistent with 
the expectation that the highest en- 
ergy galactic cosmic rays should be 
those of the highest charge. The pro- 
ton component we take to be extra- 
galactic. 

To describe the extragalactic cosmic 
rays, we assume that all sources have 
the same power law spectrum, and 
that cosmic rays lose energy in propa- 
gating to the earth by pion and e + e~ 
production from the CMBR photons, 
and by the cosmological red shift [29] . 
The sources are assumed to be uni- 
formly distributed and to evolve in 
density by (1 + z) m . Figure 9 shows 
our spectrum result with the best fit 
superimposed on it. The fitted val- 
ues m and of —7, the spectral slope 
of the spectrum at the source, are 
m = 2.6 ± 0.4 and -7 = 2.38 ± 0.05. 



10 Summary 

We have measured the flux of ul- 
trahigh energy cosmic rays from 
1.6 x 10 17 eV to over 10 20 eV. Our 
experiment detects atmospheric fluo- 
rescence light from cosmic ray show- 
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Fig. 9. E 3 times the UHECR Flux. Re- 
sults from the HiRes-I (red squares) and 
HiRes-II (black circles) detectors are 
shown. Also shown is a fit to a model de- 
scribed in the text. The la upper limits 
for two empty bins of each HiRes spec- 
tra are also shown. 

ers and performs a calorimetric mea- 
surement of cosmic ray energies. We 
perform calibrations of our detec- 
tor and measure the light-scattering 
properties of the atmosphere. The 
total systematic uncertainty in our 
spectrum measurement averages 
31%. 

In our energy range we observe two 
features in the UHECR spectrum vis- 
ible through changes in the spectral 
power law. We observe the ankle at 
3 x 10 18 eV. We also have evidence 
for a suppression at a higher energies, 
above 6 x 10 19 eV. 
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